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Abstract 

Background: Adverse remodeling of the left ventricle (LV) following nnyocardial 
infarction (Ml) leads to heart failure. Recent studies have shown that scar anisotropy 
is a determinant of cardiac function post-MI, however it remains unclear how 
changes in extracellular matrix (ECM) organization and structure contribute to 
changes in LV function. The objective of this study is to develop a model to identify 
potential mechanisms by which collagen structure and organization affect LV function 
post-MI. 

Methods: A four-region, multi-scale, cylindrical model of the post-MI LV was developed. 
The mechanical properties of the infarct region are governed by a constitutive equation 
based on the uncrimping of collagen fibers. The parameters of this constitutive 
equation include collagen orientation, angular dispersion, fiber stiffness, crimp angle, 
and density. Parametric variation of these parameters was used to elucidate the 
relationship between collagen properties and LV function. 

Results: The mathematical model of the LV revealed several factors that influenced 
cardiac function post-MI. LV function was maximized when collagen fibers were 
aligned longitudinally. Increased collagen density was also found to improve stroke 
volume for longitudinal alignments while increased fiber stiffness decreased stroke 
volume for circumferential alignments. 

Conclusions: The results suggest that cardiac function post-MI is best preserved 
through increased circumferential compliance. Further, this study identifies several 
collagen fiber-level mechanisms that could potentially regulate both infarct level 
and organ level mechanics. Improved understanding of the multi-scale relationships 
between the ECM and LV function will be beneficial in the design of new diagnostic 
and therapeutic technologies. 

Keywords: Cardiac mechanics. Myocardial infarction. Collagen fiber alignment, 
Microstructure based mechanical model. Adverse remodeling, Anisotropy 



Background 

Coronary heart disease is the leading cause of death, accounting for over 400,000 Uves 
in the United States every year [1]. Blockage of the diseased coronary arteries leads to 
myocardial infarction (MI). Post-MI, the left ventricle (LV) undergoes a complex 
remodeling process that results in the formation of a scar or infarct. The increased 
stiffness and diminished contractility of the scar reduce LV function and can lead to 
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heart failure [2]. The decrease in LV pump function has been correlated to the size of 
the infarct, but infarct size alone cannot explain patient outcome [3] . One key factor af- 
fecting LV function is the mechanical properties of the scar tissue. Increased infarct 
stiffness reduces inflation under systolic pressures, but also impairs filling of the LV 
during diastole [4], Therefore, it is imperative to elucidate the relationship between 
scar composition and mechanics in order to identify the properties that best preserve 
LV function. 

Cardiac scar tissue is composed of collagen fibers, the primary determinant of the 
mechanical properties of the scar, as well as cells and other extracellular matrix proteins 
which are considered to be a ground substance [2]. Thus, cardiac scar tissue demonstrates 
anisotropic, non-linear material behavior [5]. Naturally, the density and alignment of 
collagen fibers affect the mechanical properties of scar tissue and thus the wall stress, 
with scar tissue demonstrating a high stiffness along the collagen fiber direction and a 
lower stiffness in the cross-fiber directions [6]. Furthermore, it has been shown that 
collagen deposition during remodeling is highly controlled by the local stresses or 
strains in the tissue [7-9]. In arteries it has been shown that the angle of mean collagen 
fiber alignment adapts to optimize a tissue s load-bearing ability which leads to fibers 
that are aligned in a direction in between the directions of the two principal stresses 
[10]. In fact, recent studies by Fomovsky et al. have shown that mechanical cues are 
the primary determinant of collagen alignment [11]. Therapeutic strategies that aim to 
limit infarct expansion by reducing local stress are currently under development 
[12,13], including efforts to anisotropically reinforce the infarct region [14]. Thus, it is 
of clinical interest to determine the relationship between collagen alignment, scar tissue 
mechanics, and LV function, in order to find a possible alignment that would better 
preserve LV function. 

The objective of this study was to illustrate how collagen alignment affects LV function 
post-MI using a simple, thin-walled, cylindrical LV model with a scar region governed by 
a collagen fiber based constitutive equation. The cylindrical model is based on a concept, 
previously presented by Han et al. [15] and Oshinski et al. [16], to predict ejection fraction 
from MRI images, and the collagen fiber-based constitutive equation was based on the 
work of Grytz and Meschke [17]. This simple multi-scale model allows us to focus on 
identifying the properties of the collagen matrix that primarily determine LV function 
post-MI. 

Main text 

Methods 

Constitutive model 

The LV was modeled as a cylindrical membrane formed by folding a planar sheet 
consisting of three normal, healthy myocardial regions and a collagenous scar region 
into a cylinder (Figure 1). The mid-wall circumference of the cylinder was taken to 
be the length of the planar sheet in the circumferential direction. The healthy LV 
regions were modeled as sheets with nine layers having myocyte fiber angles ranging 
from -50 to +50 degrees relative to the circumferential direction. While reports of 
myocyte rotation in the literature are somewhat variable between species, location in the 
ventricle, and measurement technique [18-20], the distribution parameters chosen provide 
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Figure 1 Schematic of tlie four-region cylindrical model. The heart wall was modeled as blocks of 
membrane of normal contraction (R: remote) and scar (S) with roman numeral supercripts used to denote 
the three different remote regions. Labels are given to denote the circumferential length (/.), longitudinal 
height (H), and wall thickness (T) of the regions. 

V J 



a reasonable physiological approximation. The healthy region was governed by a Fung-type 
exponential strain energy density equation [21,22]: 

W = ^c{e'i-l), (1) 

with 

Q = bffE/ + b^^ {Ej + Err^) + bf^ {E/ + E/) . (2) 

Where c is a material constant b^ b^oo and bf^ are non-dimensional material constants 
and E-ff, Eca Err, Efa and £c/are the components of the Green strain tensor. The subscript/ 
denotes the fiber direction, the subscript c denotes the in-plane cross fiber direction, the 
subscript r denotes the radial (out-of-plane) direction and the subscript x denotes either 
the cross-fiber or radial direction. The Fung-type constitutive equation has been widely 
used in soft tissue biomechanics and is still commonly used for modeling myocardial 
tissue [23,24]. Since the collagen density of the healthy tissue is much lower in the 
healthy regions of the myocardium compared to the scar region, no explicit model of 
collagen is included for these regions. Cauchy stresses were determined from the strain 
energy density functions following [25], 

dW T 

M = -M/]+2[F].^.[F]^ (3) 

where is a Lagrange multiplier used to enforce incompressibility, [/] is the identity 
matrix, [F] is the deformation gradient, and [C] is the right Cauchy-Green deformation 
tensor. 

Contraction of the healthy tissue was modeled using the elastance model of active 
cardiac contraction proposed by Guccione et al. in the form of [26,27], 

(^active = (^^{Q ' (l-COSC0(4^f)) (4) 

where 4 and t are the the sarcomere length and time during a cardiac cycle, respectively. 
CO is a function of the sarcomere length and time. a° is the peak contraction stress that 
depends on cellular calcium concentration and sarcomere length [27]. The total stress 
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in the healthy regions is calculated as the sum of the active stress which is non-zero solely 
along the direction of myocyte orientation, and the passive Cauchy stress tensor. The 
individual layers of the healthy tissue were assumed to have the same planar deformation 
(i.e. the transmural shear deformation was ignored). 

The scar tissue was modeled as a thin sheet composed of a distribution of collagen 
fibers embedded within a ground matrix. The mechanical properties of the scar region 
were governed by the microstructure-based nonlinear model proposed by Grytz and 
Meschke [17,28]. Briefly, the stretch of collagen fibers was modeled as the uncrimping 
of a coiled collagen filament. This is similar to the manner in which a spring uncoils as 
it is stretched. The strain energy density function for the collagen fiber component of 
the tissue, Wcoh was calculated by integrating the first Piola-Kirchoff stress, Pcoi which is a 
function of the stiffness of the collagen filament, E^oh the ratio of the fiber diameter, Dcoi> 
to the filament diameter, dcoh and the initial crimp angle, 6o, (a measure of how tightly 
coiled the fiber is, with 0° representing a fully straightened filament and 90° representing a 
filament that has been compacted into a disc, see Figure 2) over the deformation of the 
fiber, A. 

WcoiiX) = j Pool (\ Ecol, ^ , Oo) d\ (5) 

1 \ Clcol J 

where \ depends both on the probability density function of the collagen fiber align- 
ment, /), and the right Cauchy-Green deformation tensor [C] [17,28]. The ground 
matrix was modeled as an isotropic, neo-Hookean material with a strain energy 
density 

W, = Cg{h-3) (6) 



where Cg is a material constant and Ii is the first invariant of the right Cauchy-Green 
deformation tensor. The total strain energy density, was calculated from the volume 
weighted summation of the strain energy densities of the collagen fibers, Wcoh arid the 
ground substance Wg, 

W={l-OJcol)Wg^OJcolWcol, (7) 

where Ocoi is the collagen volume fraction. For a detailed description of the calculation 
of Wcoi please refer to the original work by Grytz and Meshke [17,28]. Accordingly, 
the Cauchy stresses were determined from this strain energy density functions using 
Eq. (3). 

Boundary conditions and equilibrium equations 

For clarity, regions will be referred to as R", R"^ representing the healthy, remote 
regions, and S representing the scar region (Figure 1). In order to solve for the deformed 
volume of the LV it is necessary to establish boundary conditions and equilibrium 
equations to relate deformation to the applied loads. The base of the LV is fixed in the 
longitudinal direction, but not in the radial or circumferential directions. Further, it is 
assumed that the LV will remain cylindrical throughout deformation. Given this boundary 
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condition and the assumptions of thin walls and incompressibility, the deformation 
gradient, [F], for each region has only 3 unique components: 



Xe Kez 0 
0 0 
1 

XeXz 



0 0 



(8) 



The components Xq and represent the stretch ratios in the circumferential and 
longitudinal directions respectively and Kq^ represents the shear deformation. Since 
uniform stress distributions over each region are assumed, internal forces can be directly 
related to the Cauchy stress through the definition [25] 



= niUi ' [a] for i = 6,z, 



(9) 



where fii is the normal to the f'^ surface and Ui is the surface area of the f'^ surface. 
The internal forces must balance the forces generated by the internal pressure. Under 
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the assumptions of the Law of Laplace for thin-walled closed-end, cylindrical pressure 
vessels the magnitude of the internal forces are 

\\q\\e=PrK (lO.a) 

where qe and are the internal circumferential and longitudinal forces, P is the pressure, 
r is the deformed radius of the LV, and h is the deformed longitudinal height of the LV. 

To model the coupling of the regions and the division of load amongst adjoining 
regions, the LV wall was modeled as a 2-D system of springs (Figure 3). Based on this 
arrangement, the internal loads described in Eq. (10) are shared as 

Me = Ik/ll + Ik/'ll = Ik/ll + Ik/ll = Prh, (ll.a) 

Ik.ll = Ik/ll + Ik/ll = + Ik/ll = ^Pr'- (ii.b) 

The spring system model was also used to describe which regions share displacements. 
The deformations are described by the components of the deformation gradient tensor, 
[F], as seen in Eq, (8). In the circumferential direction, regions and R^^^ are fixed in 
parallel so that they will always have the same circumferential deformation X§, The 
same goes for regions R" and S which share circumferential deformation X^, In the 
longitudinal direction, regions R^ and R" are fixed in parallel so that they will share 
longitudinal deformation Af and shear deformation k^^- The same is true for regions 
R"^ and S that share Af, and Ko^- From Eq. (11) a full set of equilibrium equations can 
be derived which in conjunction with the constitutive equations can be solved numerically 
to determine the unknown displacements. A more complete derivation appears in the 
appendix. 
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Numerical solution methods 

The LV volume contained by each region was determined independently as: 

2 2 ' 2r 



+ (12) 



where, r is the radius of the LV, and l\ h\ and are the deformed length, height and 
thickness of each of the four regions R^^ R^^^ and S assuming that the deformed 
LV shape remains cylindrical The total volume of the ventricle was determined by 
combining the component volumes. All numerical operations were carried out using 
Matlab (Mathworks, Natick, MA). An iterative solver based on Newton's method 
with derivatives calculated by complex Taylor series expansion was used to simulate 
diastolic filling from 0 to an end diastolic pressure of 1.3 kPa (9.8 mm Hg) over 100 
pressure increments. Convergence testing showed that end-diastolic volumes calculated 
using 50 pressure increments were only 0.1% less than those calculated using 100 
increments. To simulate systole, 10 points on the end-systolic pressure volume relationship 
(ESPVR) curve were calculated using a Newtons method solution technique to solve the 
force balance equations of Eq. (11) with the added time- varying myocyte contraction forces 
and an added isovolumic constraint so that pressure could be calculated as a function of 
time. Each systolic contraction was calculated using 100 time points. From this ESPVR the 
final systolic volume was taken to be the point on the curve at 13 kPa (98 mm Hg). Stroke 
volume (SV) and ejection fraction (EF) were then calculated from the end diastolic 
and systolic volumes. 

Numerical examples 

A set of typical material properties and geometric dimensions was chosen as an example 
problem for our simulation. LV dimensions were chosen so as to roughly mimic the 
dimensions of a canine left ventricle [29]. Total LV height was chosen to be half of 
the mid-line circumference. The thickness of the infarct was chosen to be half the 
thickness of the healthy region as well as half of the total circumference and half of 
the total height. In this configuration the scar represents 25% of the total myocardial 
area representing a moderately sized scar. The final dimensions were chosen so that 
the LV volume in the undeformed state matched that reported by Fomovsky et al. 
[24]. The material constants of the Fung equation were taken from Fomovsky et al. 
[24], with incompressibility assumed. The properties governing myocyte contraction were 
the same as those reported by Guccione, and an initial sarcomere length of 1.65 (im was 
chosen as this resulted in a baseline ESPVR that roughly matched the baseline case 
reported by Fomovsky et al. [24]. The collagen fibers of the scar region were aligned 
following a Von Mises distribution wrapped between -90° and +90°: 

exp(y* cos(2(p)) 
nl.(y) 

where y is z. fiber alignment factor such that high values of y correspond to highly 
aligned fibers, and 4 is the modified Bessel function of order zero. The properties of 
the scar region were determined by assuming a collagen fiber density of 40% and a 
circumferential mean fiber alignment and fitting the remaining parameters (crimp 
angle, young s modulus, neo-Hookean constant, ratio of helix diameter to filament 
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diameter, and fiber alignment parameter) to the stress strain relationship reported by 
Gupta et al [5] for a 6 week old infarct in a porcine model A circumferential alignment 
was chosen as this was the axis of the largest stress under equibiaxial stretch in the 
experimental data. The values for these parameters are listed in Table 1. Least-square 
fitting was performed using a gradient based method and the value was 0.996 for 
the circumferential stress and 0.992 for the longitudinal stress (Figure 4). Using the 
properties described in Table 1 as a reference point, parametric studies were run to deter- 
mine the effects of varying the mean collagen alignment angle, collagen density, collagen 
fiber alignment, collagen fiber stiffness, collagen crimp angle, scar size, and scar thickness. 
Additionally, a baseline case with 4 healthy regions was run for comparison. 

Results 

The SV and EF of the baseline (pre-MI) case were 22.0 mL and 0.486 respectively. In 
order to investigate the effect of varying the mean fiber angle, an example infarcted LV 
was run 19 times with the mean fiber angle ranging from 0° (circumferential alignment) 
to 90° (longitudinal alignment). Simulation results showed that SV and EF increased 
with fiber angle and the maximum SV and EF were achieved when the fibers were 
aligned at 90° (Figure 5A-B). SV and EF were reduced compared to baseline through 
lowered diastolic filling and a right shift of the ESPVR (Figure 5C). Longitudinally 
aligned fibers improved diastolic filling when compared to circumferentially aligned fibers. 
This extra diastolic deformation also allowed myocytes in the healthy regions to generate 
more contractile force which slightly reduced end-systolic volume as seen in the ESPVR s 
in Figure 5C. At angles intermediate to 0° and 90°, SV and EF rise steeply from around 20° 
to 60° before plateauing. 

Table 1 Parameters of the example infarcted LV model 



Geometry 

Total undeformed LV length 10 cm 

Total undeformed LV height 5 cm 

Undeformed thickness of healthy regions 0.8 cm 

Undeformed thickness of scar region 0.4 cm 
Mechanical properties of healthy regions 

Fung constant (c) 880 Pa 

Fung constant (bcc) 18.5 

Fung constant (bxx) 3.58 

Fung constant (bxc) 1.63 

Number of layers 9 

Myocyte fiber angle range -50 to 50 
Mechanical properties of scar region 

Collagen crimp angle (Go) 25.5° 

Ratio of collagen helix diameter (Dcoi) to filament diameter (dcoi) 2.39 

Young's modulus of collagen (Ecoi) 1.16 MPa 

Neo-Hookean constant for ground material (Cg) 5.72 kPa 

Fiber alignment factor ( Y ) 0.771 

Volumetric collagen fiber fraction (cjOcoi) 0.4 
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Parameter Fitting: Circumferential Parameter Fitting: Longitudinal 




Circumferential Stretch Ratio Longitudinal Stretch Ratio 

Figure 4 Fitting of material properties to experimental data. The properties of the scar tissue were fit 
to experimental data obtained by Gupta, et al. 1994 for a 6-week old porcine infarct. (A) The fitting to the 
circumferential stress resulted in an value of 0.996. (B) The fitting to the longitudinal stress resulted in an 
value of 0.994. 



The deformations of the scar and remote regions (S and R^), as measured by the 
components of the deformation gradient (see Eq. (8)), were also examined as the mean 
collagen angle was varied (Figure 6). When collagen aligns towards the longitudinal 
direction (90°), the circumferential deformation of the scar region increased while the 
longitudinal deformation at end-diastole slightly decreased (Figure 6A). Longitudinal fiber 
alignment also reduced the longitudinal deformation of the region while having only a 
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Figure 5 Cardiac function of an example infarcted LV. Changes in EF (A) and SV (B) with the mean 
alignment angle of the collagen in scar tissue. Dotted lines are the baseline values corresponding to the 
healthy LV. Both EF and SV were maximized when collagen fibers were aligned longitudinally. The pressure 
volume relationships (C) show that SV is reduced through both reduced diastolic filling and a right shift in 
the ESPVR. Longitudinal alignment of fibers reduces both of these effects. 
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Figure 6 Deformation and stress. Diastolic deformation of the scar (A) and remote region R (B) and 
systolic deformations of the scar (C) and R' regions (D) plotted as functions of mean collagen alignment 
angle. The normal components of the deformation gradient (circumferential and longitudinal) are shown 
on the left vertical axis and the shear component is shown on the right axis. 



slight effect on the circumferential deformation at end-diastole (Figure 6B). At end-systole, 
longitudinal fiber alignment greatly reduced longitudinal deformation in the scar region 
while only slightly increasing the circumferential deformation (Figure 6C). Longitudinal 
fiber alignment reduced longitudinal deformation in the remote region but also slightly 
increased the circumferential deformation at end-systole (Figure 6D). 

Several parametric variations in material properties and geometry were simulated and 
are presented as contour plots (Figure 7). Increasing the size of the scar (as measured 
by the percentage of the total LV area greatly reduced SV (Figure 7A). The effect of 
increased longitudinal alignment was greatest for scars that were about 20-30% of the 
total LV area. Reducing the thickness of the scar tended to improve SV by allowing for 
greater diastolic deformation (Figure 7B). However, decreasing scar thickness below 
20% of the thickness of the healthy regions resulted in decreased SV. Increasing the 
collagen fiber density improved SV for longitudinally aligned fibers, due primarily to 
increased diastolic deformation in the cross-fiber direction (Figure 7C). For circumferential 
alignments, increased collagen density only slightly reduced stroke volume. Increasing the 
neo-Hookean constant governing the ground substance of the scar region reduced the SV 
(Figure 7D). This effect was larger for longitudinally aligned fibers. When the Youngs 
modulus of the filaments was low, SV was correspondingly low (Figure 7E). Optimal 
stroke volumes were found when fibers were aligned longitudinally and the Young s 
modulus was about 1 MPa. Increasing the Youngs modulus to 10 MPa tended to reduce 
SV, however increasing the modulus further to about 100 MPa increased the SV. Varying 
the initial crimp angle revealed that optimal SV occurred at around 20 degrees (Figure 7F). 
Reducing the alignment factor y, a change that represents a more disperse distribution of 
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Figure 7 Parametric study of SV. Contour plots of stroke volume are presented with collagen alignment 
angle on the x-axis and one of the following parameters of interest on the y-axis. (A) Infarct size, presented 
as percent area of total LV area. (B) Wall thickness, presented as the ratio of the thickness of the infarct 
region to the thickness of the remote region in percent. (C) Collagen density presented as the percent 
of scar volume. (D) The log of the neo-hookean constant for the ground material of the infarct in log 
(Pa). (E) The log of the Young's modulus of the collagen filament in log(Pa). (F) The initial crimp angle 
of the collagen filament in degree. (G) The log of the fiber alignment factor y. (H) To aid in the interpretation of 
the contour plots, a plot of stroke volume vs. mean collagen alignment angle is presented for three levels of 
the fiber alignment factor y. 



fibers, increased the SV w^hen the mean angle was near 0° but reduced the SV w^hen the 
angle was near 90° (Figure 7G and H). 



Discussion 

We developed a cylindrical LV model that links the wall stress and pumping function 
of the dysfunctional LV to the density and alignment of collagen fibers in the scar tissue. 
Using this model, we examined the effect of collagen alignment and density on LV 
function post-ML Our results demonstrate that the orientation of collagen fibers in 
the scar region of the myocardium has a large influence on the performance of the LV. 
It was found that longitudinal fiber alignment limits scar deformation during systole 
and maximizes LV performance. 

The current model clearly demonstrates how the mechanics of the healthy contractile 
tissue are coupled to the properties of the collagen network in the scar region. The 
conclusion of longitudinal aligned fibers maximizing the LV function agrees with reports 
in the literature [24]. This study also ellucidates the mechanism behind this improvement. 
The results suggest that rather than increased longitudinal stiffness improving function it 
is actually increased circumferential compliance that has the largest benefit. Increased 
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circumferential compliance allows for the greatest diastolic filling, which in turn allows 
for greater contractile ability through the Frank-Starling mechanism. It is worth noting 
that an infarct that is too compliant may eventually lead to reduced systolic function as it 
inflates under systolic pressure. Examining the ESPVRs in Figure 5C, it is seen that the 
longitudinal and circumferential alignment curves cross once the pressure exceeds about 
16 kPa indicating that the increased circumferential compliance is leading to increased 
infarct bulging under high pressures. 

While the parametric studies demonstrated that smaller scars preserve LV function 
they also revealed that for certain scar sizes, thinner scars increase LV function, a far less 
intuitive conclusion. This is due to the fact that thinner scars allow for more diastolic 
filling. There is a limit to this though, as scars with an infarct thickness to remote 
thickness ratio of less than 20% had a reduced SV due to increased systolic bulging. 
We do not endorse infarct thinning as an advantageous phenomenon as it increases 
the risk of cardiac rupture and infarct expansion. 

The finding that increased collagen density has a positive effect on SV for longitudinally 
aligned fibers but very little effect for circumferential alignments is explained by the 
fact that as more collagen fibers replace the ground substance, the tissue becomes 
more compliant in the cross-fiber direction but stiffer in the direction of the fibers, 
which leads to increased diastolic filling. The importance of increased compliance is 
further seen in the findings that increasing the neo-Hookean constant decreased SV. 
One important takeaway from the parametric studies is that there appears to be an 
ideal stiffness at which diastolic filling is increased but systolic contraction is not 
greatly impaired. This is seen in the results for the infract thickness. Young s modulus, 
and initial crimp angle. It is also interesting that the fitted values we found for the 
Young s modulus and the crimp angle of porcine scar tissue appear to be near these ideal 
stiffness regions. The finding that increased fiber dispersion is beneficial for circumferen- 
tially aligned fibers but detrimental for longitudinally aligned fibers further highlights the 
importance of scar anisotropy in preserving ventricular function. 

This work was strongly motivated by the work from the Holmes group [24]. They 
studied the effect of scar anisotropy on LV function using a finite element model with 
mechanical properties of both the healthy and infarct regions described by a Fung-type 
constitutive equation. By comparing isotropic scars, longitudinally stiff scars, and circum- 
ferentially stiff scars, they found that SV would be maximized when the scar tissue was 
stiffer in the longitudinal direction than in the circumferential direction. Fomovsky s work 
identified the direct role that infarct stiffness plays in controlling LV mechanics so as to 
inform the design of infarct restraint or stiffening therapies. Our current study developed 
a simple model to relate cardiac function to the features of the scar extracellular matrix 
rather than simply to stiffness through the use of a microstructure based constitutive 
model. The agreements that we see between our results and the results of the more 
complex finite element model presented by Fomovsky, in terms of increased longitu- 
dinal stiffness leading to increased SV give us confidence in our findings. A more 
recent study from the Holmes group has demonstrated that regional mechanics control 
collagen alignment [11]. In complement, we have demonstrated here that collagen align- 
ment affects the scar mechanical properties and thus LV function. The model we have 
presented could be a simple, easy to use, analytical tool for studying how initial collagen 
alignment controls remodeling and infarct expansion through mechanical feedback. 
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While this model has offered much insight into the importance of collagen alignment 
in determining heart function, there are several limitations to the study. This model 
assumes a thin walled geometry which ignores the effects of transmural variation in 
stresses. While the choice of the thin-walled assumption may create additional error in 
the calculated stresses, it was deemed acceptable since the focus of the study was on 
identifying potential mechanisms rather than conducting detailed local stress analysis. 
Based on the fact that infarct scars are thinner than normal LV wall, we chose to model 
the infarct as a single layer with a disperse fiber alignment. In reality the infarct is 
composed of several layers of collagen fibers with a mean fiber orientation that varies 
through the wall thickness similar to the helical pattern of myocyte alignment in normal 
myocardium. This simplification in itself should have little effect on the results due to 
the thin wall assumption of the model. Further, the choice of a single fiber distribution 
simplified the fitting of the scar properties to experimental data. This model also 
neglects the effects of ventricular torsion, which plays an important role in systolic 
ejection [30]. 

In addition, the assumption that adjacent regions share deformations along the shared 
edges allowed for a simple method of enforcing material compatibility. On the downside, 
this means that the stress field is non-continuous. This assumption was made to capture 
the phenomenon of infarct coupling whereby increasing the stiffness of the infarct reduces 
the deformation of the neighboring regions. In reality this effect is much smaller in 
regions that are far from the infarct region than in the small border zones that surround 
the infarct with properties that transition from scar to normal, making it possible that the 
magnitude of the fiber orientation effects are slightly overstated by the model. Given that 
the properties in the normal regions should be independent of the infarct properties 
we believe that the mechanisms through which scar architecture influence SV identified 
by this model are correct. In whole we believe that the model is not oversimplified as it 
matches behavior reported both by Fomovslcy, (improved SV with increased longitudinal 
stiffness) and by Bogen (reduced stiffness allowing for greater diastolic filling and contrac- 
tion for small infarcts) [4,24]. 

The ground-up modeling approach we have chosen wiU allow us to refine the model 
in future work. While the collagen network is the primary determinant of infarct 
mechanics in the fibrotic scar, there are other structural proteins such as elastin and 
ECM modifications such as cross-linking that also play a role in LV mechanics [2]. 
Refinements of the constitutive model to include the effects of these ECM features 
would be an interesting direction for future work. In addition, the ischemic LV is 
highly inhomogeneous and further effort to model the border zones as well as the 
transmural variations in infarct properties would add to the physiological relevance of the 
model. Another future direction could examine the mechanics of the LV during the ische- 
mic or inflammatory phases occurring immediately following MI. During these phases 
mechanical properties of the infarct are controlled less by collagen and more by the 
turnover of existing ECM and intracellular structural proteins such as titin as weU as 
edema. The development of the collagen network and in turn the mechanical properties 
of the scar following this initial period are highly controlled by the expression of growth 
factors such as TGF-p and proteases such as MMPs [3,31]. Coupling the mechanical 
model to a dynamic model of protein regulation is one potential direction for future 
research. Another dynamic process that could be modeled in future work is the change 
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in the orientation of the fibers that occur under load, which can dynamically change 
the mechanical properties of the tissue [32,33]. 

The post-MI LV is a complex and dynamic physiological environment and much 
work is ongoing to create highly accurate models. The goal of our paper was not to 
improve upon the accuracy of the existing numerical models, but to instead clarify the 
role that collagen plays in controlling infarct mechanics. Compared to some of the 
complex finite element models found in the literature, the simplicity of this model is 
actually an advantage as it focuses attention to the most important variables. This 
makes it a quick tool for identifying the key structural determinants of LV function 
post-MI. For example, the current simplicity of the model may make it a good choice 
for investigating many of the temporal changes that occur post-MI. 

In summary, our model results showed that EF and SV are both optimized when fibers 
are aligned longitudinally, which further supports findings in the literature. The finding of 
an ideal fiber angle may aid in the development of diagnostic indices and treatment strat- 
egies based on determining fiber angle (perhaps through angle-sensitive MRI [34]) and 
directing it towards an optimal orientation through controlled remodeling. Since the 
model is based on the mechanical properties of actual infarct tissue it will also serve as an 
excellent starting point for creating multi-scale models that explore the interactions of gen- 
ome, proteome, cell, tissue, and organ level changes in the post-MI remodeling process. 

Conclusions 

A simple model of the post-MI LV was developed to demonstrate that a longitudinal align- 
ment of collagen fibers in the scar region maximizes cardiac function. This improvement is 
the result of increased compliance in the circumferential direction which allows for greater 
diastolic filling and greater systolic contraction through the Frank-Starling mechanism. 

Appendix 

1. Derivation of equilibrium equations 

Many of the basic continuum mechanics equations used below can be found in Humphreys 
Cardiovascular Solid Mechanics text [25]. As described in the paper, there are six unknown 
deformation components. In order to solve for these six deformation components it is 
necessary to derive six independent equilibrium equations to relate stress to deformation. 
The internal forces in each region that balance the intraventricular pressure can be calcu- 
lated from the Cauchy stress given in Eq. (9). If a matrix [a] is defined with rows njui and 
a matrix [^] is defined with rows qi then Eq. (9) can be written as 

m = [a]-[a]. (Al) 

Similarly a matrix [A] can be defined with rows NiAi representing the areas in the 
undeformed configuration and a matrix [E] with columns Qi representing the forces in 
the undeformed configuration. Both [S] and [A] can be related to [^] and [a] respectively 
through the deformation gradient. 

m = mFf, (A2) 

[a] = [A]-[F]-\ (A3) 

In the undeformed configuration, it is assumed that the surface normals Ni are aligned 
with the circumferential, longitudinal, and radial axes and that the internal forces in the 
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undeformed configuration Qi are normal to the undeformed surface normals Ni since 
pressure acts normal to surfaces. Thus, 



(A4) 





'Qee 
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0" 
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Qzz 
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0_ 
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~HT 
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LT 
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LH 



(A5) 



where, H is the undeformed, longitudinal height of the region, L is the undeformed, circum- 
ferential length, and T is the thickness of the region. Now solving Eq. (A2) in matrix form. 



Qee 0 0 
0 Q.. 0 

0 0 0 



0 0 

Kez 0 
0 0 



1 

XbXz 



0 0 
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0 0 0 



^ee 0 0 
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0 0 0 



(A6) 



it is seen that there are only 3 non-zero internal force components of [^]. Eq. (A3) can 
also be solved in matrix form for [a] 
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(A7) 



Substituting Eq. (A7) into Eq. (Al) and assuming negligible radial stresses yields 
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(A8) 

Relating the result of Eq. (A8) with the components of [^] as given by Eq. (A6) provides 
a system of 4 equations 



(^66 -(^zd 



Kez\ HT 



0 = ( (Tdz- (^z. 



LT 

Xz 

LT 



Xz J Xe 
Kez\ HT 



Xz J Xq 



(A9.a) 
(A9.b) 
(A9.c) 
(A9.d) 



From Eq. (A9.b), it is seen that 



KQz 



(AlO) 
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Combining this fact with Eq. (A9x) and Eq. (A9.d) yields the following relationship 
between force and deformation 

^ = ^^. (All) 

Now that relationships between force, stress, and deformation have been established 
for individual regions it is necessary to examine the four region model with forces and 
deformations being shared as described in Figure 3. In order to determine the magnitudes 
of the internal force vectors described in Eq. (11) it is necessary to examine the components 
of these vectors again. 
From Eq. (A6) it is seen that 
^6 = We^ (A12.a) 
^z = We + (Izzh, (A12.b) 
where is and 4 are unit vectors. From Eq. (A12.a) it is seen that 

he\\=%e^ {A13) 
which allows for Eq. (11. a) to be rewritten as 

^do = (lee + (lee^^ = ^ee" + ^ee^ = P^h. (A14) 
Substituting Eq. (A9.a) for each region into Eq. (A14) yields 2 equilibrium equations 

cree' - oj ^ + f <^ee''' - crj" ^ ) '-^^ = P^K (A15.a) 





f^iiij^iii 




1 J /// 

A0 






' A/y 


' A/ - 



^z^ J ^e 

e" - oJ' ^) + ioee' - a,e' ^ ) = Prh. (AlS.b) 

Az J Ae \ Az J Aq 

The vector can be written in polar notation as 

^z = lkzl|sin(c|))«0 + ||^J|cos(c|))4. (A16) 
where 

(|) = atan(^^^. (A17) 

Substituting Eq. (All) into Eq. (A17) yields 
(|) = atan(^^^. (A18) 

Equating Eq. (A16) with Eq. (A12.b) gives 
^ze = lkzl|sin(c|)), (A19.a) 
q^^ = \\qz\\cos{cp), (A19.b) 

Solving Eq. (ll.b) for the magnitude and equating with Eq. (A19) yields 

II 11 ^z6> ^ ^z/ I ^z/^ ^ ^z/^^ I ^ze^ ^ p 2 (A1(^^^ 

sin(c|)) sin((|)^) ^ sin(c|)^^) sin(c|)^^^) ^ sin(c|)^^) ' ^ ^ ^ 



/ ^11 ^ III ^ s 

nPr^. (A20.b) 



^zz ^zz , ^zz ^zz , ^z _n,.2 



cos((|)) cos((|)^) cos((|)^^) cos((|)^^^) cos((|)^^) 
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Substitution of Eq. (A9x) for each region into Eq. (A20.a) gives 2 more equilibrium 
equations 

TT-TTTT + "o" ,11 ■ t = ' ^^^^^ 
Az sm((|)^) Az sm((|)^^) 



jlllq-III T^T^ 
III ^ ^ 1^5^^ 



+ ^^-^ , c . , = T^Pr". (A21.b) 



A/'sin((|)^^^) ' A/sin((|)^) 



Substitution of Eq. (A9.d) for each region into Eq. (A20.b) gives the last 2 equilibrium 
equations 

a J —r h (Tzz —Tt = T^Pri (A22.a) 

A/cos((|)0 Xj'cos{^'') ' ^ ^ 

(Tzz"^ —m + (^zz^ = ^Pr^ • ( A22.b) 

Xj"cos{(p"') A/cos((|)^) ^ ^ 

Collecting Eq. (A15), Eq. (A21), and Eq. (A22) and noting the shared lengths, heights, 
and deformations as shown in Figure 1 produces the final set of equilibrium equations. 

/ , rKez''\H^T' ( jjjKez^\H^T'" 
(Tee - (Tze -f^ — ^ + oee"' - Ozo"' -fs = (^^S.a) 

V Az / A0 V Xz / A0 

(^66 - (^zd . 5 ^ \ (^dd - (^zd , s = (A23.bj 

Az / A^ V Az / Aq 

jR jl jS jll 
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Az Az 
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Az Xz 
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Ozz^^ — ^ + Ozz^ — ^ = TTPr^cos ((|)^) , ( A23.f ) 

Az Xz 

where 



R 

c|)^ = atan 



c|)^ = atan(^^^. (A24.b) 

This set of 6 independent equilibrium equations can be used in conjunction with the 
constitutive equations to solve for the 6 independent deformation components, A^,A;^, 

KXAz^^^Az' 
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